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ABSTRACT 



Small-scale magnetic fields are thought to play an important role in the heating of the outer solar atmosphere. By taking advantage 
of the unprecedented high-spatial and temporal cadence of IMaX, the filter vector polarimeter on board the Sunrise balloon-borne 
observatory, we study the transversal and longitudinal velocity oscillations in small magnetic elements. The results of this analysis are 
then compared to MHD simulations, showing excellent agreement. We found buffeting-induced transverse oscillations with velocity 
amplitudes of the order of 1 - 2 km/s, to be common along with longitudinal oscillations with amplitudes ~ 0.4 km/s. Moreover, 
we also found an interaction between transverse oscillations and longitudinal velocity oscillations, showing a ±90° phase lag at the 
frequency at which they exhibit the maximum coherence in the power spectrum. Our results are consistent with the theoretical picture 
in which MHD longitudinal waves are excited inside small magnetic elements as a response of the flux tube to the forcing action of 
the granular flows. 

Key words. Sun: photosphere. Sun: oscillations. Sun: helioseismology 
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1. Introduction 

Magnetic fields play a major role in the dynamics and the 
energetics of the solar atmos phere. They are pre sent over a 
wide range of spatial scales ("So lanki et al.l \200 6). from very 
small elements (Lagg et al. 2010), at or below the best current 
spatial resolution achieved by modem solar telescopes (^ 100 
km), up to large sunspots with diameters of order of 70000 km. 
Interestingly, magnetic fields, organized in structures similar 
to flux tubes or sheets, can act as a wave-guide for waves and 
perturbations, through differen t layers of the solar atm osphere 
JNarain & Ulmschneided[l996t iDe Pontieu eTall |2004 . It has 
been estimated that concentrated small-sc ale magnetic field s 
cover roughly 1% of the quiet solar surface (iBonet et al.ll2012h . 
Much work has been done, both theoretically and through 
numerical simulations, to investigate and to model the 
wave propagation in sl ender magnetic fl ux tu b es (e.g. 

"1911 'Roberts' 
1998yHasan et al. 



Roberts & We bbI 119781; 'Edwin & Robertj 
1983; MusielaketaLnl989; Steineretal 



2003t jMusielak & Ul mschneided l2003at 

iFedun et al.. .201 U to name a few). In contrast, com- 



Khomenko et al.l 



paratively few observational studies exist about waves in the 
abundant small-scale magnetic elements in the solar photosphere 
(Volkmeret al. 1995; Jess et al. 2009; Martinez Gonzalez et al. 
l201lt IJess et al.l l2012allb ); although oscillations and waves in 
faculae and network regions at lower spatial r esolution have 
been studied more (iKhomenko et al.l l2008al ; ICenteno et alj 
120091) . 

The small-scale magnetic features are expected to harbour 
a rich variety of waves excited largely by the buffeting they 
experience at the hands of the granulation and partly their 



interaction with the ubiquitous p-modes. iHasan et aP (l2003h 
argued that horizontal motion of magnetic elements in the 
photosphere can generate enough wave energy to heat the 
magnetized chromosphere. They found, based on numerical 
modelling, that granular buffeting excites kink waves and, 
through mode coupling, lo ngitudinal waves. This was confirmed 
by iMusielak & UlmschneideJ (l2003bi) . They in fact argued 
that although transverse tube waves do not generate observable 
Doppler signals, when observing vertical flux tubes at disk 
centre, they excite forced and free longitudinal oscillations 
through a non-linear coupling. These should be observable if 
data with sufficiently high spatial resolution and sensitivity are 
available. 

From the observational point-of-view, IVoIkmer et al.l (Il995h . 
using high spatial and temporal resolution spectropolarimetric 
data, detected short-period longitudinal waves (P ^ 100 s) 
in small magnetic elements in the solar photosphere and 
estimated the energy flux they carried to be sufficient for the 
heating of the bright structures ob s erved i n the chromo spheric 
network. Mart inez Gonzalez et al.l (1201 ll) have found, using 
SUNRISE/IMaX data, magnetic flux density oscillations in 
internetwork magnetic elements, which they interpreted to be 
due to granular forcing. 

In this work, we study the longitudinal and transversal oscilla- 
tions of small magnetic elements in the solar photosphere and 
their power spectra. To do this, we exploit the unprecedented 
combination of high spatial and temporal stability provided by 
SUNRISE/IMaX, together with the large number of observed 
magnetic features, to infer statistically significant information 
on oscillations of small-scale flux-tubes and the interaction 
between longitudinal waves and transverse kink waves. A 
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particular advantage of this approach is the prove n ability of 
SUNRISE to resolve small magnetic elements ('Lag g et all 
I2OIO), a nd to detect some of the intern al structure of network 
features (JMartinez Gonzalez et al.ll2012l) . 

This paper is organized as follows. In section 2 we outline the 
data set used in this work, as well as the MHD simulation used 
to check the validity of our observational results. In section 3, 
we describe the method used for the analysis of the oscillations 
which is based upon FFT analysis and wavelet analysis. In 
section 4, we present the results. Starting from the FFT power 
spectra obtained from a limited number of magnetic features, 
we then apply wavelet analysis to a single magnetic element 
as a case study and then extend this analysis to the full sample 
of magnetic features at our disposal, to obtain statistically 
significant phase lag information between the longitudinal and 
transversal velocity. Section 5 is devoted to the discussion of the 
results, while our conclusions are given in section 6. 
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Fig. 2. Histograms of size of the analysed magnetic elements in the 
IMaX data (continuous line) and the simulation (dashed line). The size 
of each magnetic element is defined as the equivalent diameter associ- 
ated to a circle with the same area as the measured one. 



2. Data sets employed: SUNRISE/IMaX 
observations and MHD simulations 

The data set used in this work consists of a 2D spectropolari- 
metric time series with a length of approximately 32 min and 
a cadence of 33 s, acquired by the Im aging Magnetograph ex- 
periment (IMaX; Martinez Pillet et al. 2011) on board the SUN- 
RISE balloon-borne mission (Ba rthol et al. 201 0) in the Fe I 
525.02 nm spectral line. The data were taken on 2009 June 9 
and encompass a quiet Sun region of approximately 40 x 40 
arcsec close to disk center In addition to c ommon calibrations , 
the data were phase-diver sity reconstructed (Sol anki et al.ll2010l: 
iMartinez Pillet et aDl201 1.) . resulting in a spatial resolution of 
0.15 - 0.18 arc sec. The calibrati on procedure of IMaX data is 
summarized bv lRoth et al.l (l2010l) . although we used spectropo- 
larimetric inversions instead of a gaussian fit to estimate the 
Doppler velocity in small magnetic elements. There is an intrin- 
sic advantage in doing this. In magnetic elements, the gaussian 
fit can often fail as the Stokes-I signal is reduced. The inver- 
sions, however, also include Stokes-V profiles, which are in turn 
very strong. The zero crossing wavelength point of Stokes-V is 
a velocity indicator available to the inversion code but not to the 
gaussian fits. For more details on the spectropolarimetiic inver- 
sions we refer to iGupIielmino et al.l (J2012). In Fig. \T\ we plot 
the continuum intensity map (panel a) and the longitudinal mag- 
netic field (pan el b) estimated by means of SIR spectr opolarimet- 
ric inversions (iRuiz Cobo & del Toro Iniestal 1 19921) . assuming 
a single-component atmosphere with height-independent mag- 
netic vector 

The MHD simulations analysed here were carried out using the 
MURaM code which solves the compressible MHD equations 
with an energy equation including radiation transfer in a non- 
grey approximation, and an equation of state which includes the 
effects of partial ionization (see Vogler et al. 2005, for a full de- 
scription of the equations and the numerical detail s). The code 
has been used in for both quie t-Sun studies (e.g. iKeller et al] 
l2004t IVogler & Schussledl2007l) as well in studies o f magnetic 
structures (for example pores, ICameron et al.' (TOOT), sunspots 
[kempel et al. (2009), and active regions Cheung et al. (2007) 
amongst others). Here we use a simulation domain which is 24 
Mm in both horizontal directions and 1 .4 Mm in the vertical di- 
rection. The grid spacing is 20.8 km in both horizontal direc- 
tions, and 14 km in the vertical. The top boundary is open to 
flows, and the magnetic field there is forced to be vertical. In this 
study we only use the vertical magnetic field and vertical veloc- 
ity from the tross - 1 surface (at a height of about 800 km above 



the bottom of the box), which was stored every 5 to 7 seconds. 
The total duration is 58 minutes in total. 



3. Methods 

3.1. Features tracking 

With the aim of studying the properties of the oscillations in 
small magnetic structures, we tracked the magnetic features in 
both the simulation and the IMaX data using the YAFTA code 
(Yet Another Feature Tracking Al gorithm); a labelling flux- 
ranked uphill gradie nt algorithm (IWelsch & Longcopel l2003t 
iDeForest et al.ll2007h . The algorithm tracks and labels groups of 
pixels lying on the same 'hill' in circular polarization maps, in 
the case of IMaX data, and vertical magnetic field maps in the 
case of simulations. For this purpose we set up two thresholds 
to avoid spurious detections in the final results. The first detec- 
tion threshold is set at 2cr to both Stokes V maps and B. maps, 
where for the observations the cr value corresponds to the noise 
at the IMaX continuum wavelength, whereas for the simulations 
it coiTesponds to the standard deviation of the magnetic signal 
over the FoV. It is worth to mention here that the choice of cr 
may aff'ect the final results of our investigation. For this reason 
we also did a sensitivity test using different thresholds (see Sect. 
4.3). The second threshold acts spatially, allowing the detection 
of only those magnetic features whose area is larger than 9 pixels 
(slightly larger than the spatial resolution achieved) in the case 
of IMaX data, and 49 pixels for the MHD simulation. These two 
values correspond to a linear spatial scale of 150 km and 140 
km respectively. This thresholds have been chosen in such a way 
that the distributions of size of the magnetic elements, in the sim- 
ulation and in the IMaX data, are comparable (see Fig. |2]l. The 
size of each magnetic element is computed as the equivalent di- 
ameter of the circle with the same area as the measured one. The 
position of each tracked magnetic element is obtained from the 
center of gravity, and from t hat the velocity is thus es timated. 
As extensively discussed in iJafarzadeh et al.l (l20I3h . the quan- 
tity |Bcoi(y)|, with 7 being the magnetic field inclination, is re- 
trieved more reliably than B alone, since the latter can be af- 
fected by the noise in Stokes Q and U. For this reason the use 
of Stokes V maps in the case of IMaX observations is pre- 
ferred. We are aware, however, that using B, for the simula- 
tions and Stokes V for the observations does introduce a cer- 
tain bias. Thus, we expect that we choose fewer strong-field fea- 
tures in the observations, since these generally show a significant 
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(d) IMaX features 

Fig. 1. (a) IMaX continuum image, (b) IMaX longitudinal field obtained from SIR inversions, (c) MHD simulation: longitudinal component of the 
magnetic field, (d) IMaX features tracking example. The region between the two concentric circles, in the lower right part of the image, illustrates 
the size of the aperture placed around each magnetic element to deteiTnine the properties of the p-mode oscillations near that magnetic element. 
Each coloured contour represents a label given to each identified magnetic feature (see main text). The inset displays a blow-up of a part of the 
field of view, including the left part of the large network patch near the bottom of the full image. 



line-weakening due to the high temperature (e.g. ISolankil [19861; 
iLagg et al.ll201 0). 

The visual inspection of the results of the tracking revealed a 
number of short-lived magnetic features which were observable 
(traceable) for only a few time steps. Therefore, an additional 
temporal threshold was added to exclude them from the analy- 
sis, since we are mainly interested in improving the frequency 
resolution in the spectral domain. All the features lasting for less 
than 150 s in both the simulation and the observational data set 
were rejected. In Fig.[T](panel d) we mark the features fulfilling 
the criteria in a SUNRISE/IMaX snapshot by coloured contours. 
Each color represents a label given to the identified features. 
The inset highlights that the algorithm tends to break up larger 
magnetic features into smaller ones, according to the number of 
maxima of the signal present there. 



3.2. Analysis of the oscillations 

For each magnetic element tracked by the algorithm, we esti- 
mated the vertical velocity (v^) within its area, by averaging over 
values from each individual pixel, and the horizontal velocity 
(v;,) obtained by following its position. The first one corresponds 
closely to LoS velocity in the case of IMaX data, since the FoV 
is very close to disk center The Doppler velocity was taken from 
the SIR inversions. In addition, we also estimated the contribu- 
tion of the local non-magnetic oscillatory field (vlnm) by con- 
sidering the surroundings of each magnetic element, to further 
check the results. To this end, we surrounded each magnetic el- 
ement with an aperture like the one shown in Fig. [T](lower right 



part of panel d) between the concentric circles with inner and 
outer radii of 380 km and 950 km, respectively. The ambient os- 
cillatory signal vlnm is estimated by taking the average of the 
velocity of the non-magnetic pixels (below 10 G) within the de- 
fined aperture. This prevents the contamination by possible mag- 
netic field effects. 

Our analysis consists of three steps. At first, we selected suit- 
able magnetic features, on the basis of their lifetimes, to estimate 
the average power spectral density for both longitudinal and 
transversal velocity oscillations. We chose twenty cases for both 
simulations and observations with similar long lifetimes. We re- 
stricted ourselves to a few examples, picking particularly long- 
lived magnetic features, that can harbour also lower frequency 
waves, unlike the short-lived majority of the magnetic features. 
For each magnetic feat ure we estimated the peri odogram us- 
ing the Welch method ( Blackman & Tukevlll958h . Each peri- 
odogram is then interpolated over the same frequency grid, to 
take into account for the different lengths of the time series used, 
although the lengths of the time series chosen was very simi- 
lar (on average ~ 12-13 min). We then estimated the average 
power spectral density for both longitudinal and transversal os- 
cillations. In this case the simulation was resampled to the IMaX 
cadence (33 s). 

In the second step of our analysis, we took the longest lived 
magnetic feature, among the twenty chosen before, as a case 
study. The amplitude of the oscillations of v- and v/,, as well 
as their coherence and phase were studied using wavelets. 
Since the time series associated to the transverse and longi- 
tudinal velocity are expected to be non-stationary, the wavelet 
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Fig. 3. (a) Power spectra of v/, obtained from IMaX (solid line) and simulations (dashed line). The dot-dashed line represents the average power 
spectrum of 500 simulated random walks, (b) Power spectra of v- obtained from IMaX and simulations. Each spectrum is the average of spectra 
of twenty individual features. The vertical continuous and dashed lines indicate the lowest frequency due to the total length of the time series for 
the vertical velocity and the horizontal velocity, respectively, (c) Amplitudes of v;,. (d) Amplitudes of v-. The vertical continuous and dashed lines 
indicate the mean value for IMaX and the simulation respectively. 



analysis presents many benefits over the Fourier-based anal- 
ysis, in particular , with respect to the estimate of the phase 
(Blo omfield et alJ 2004). O ur analysis was performed with the 
standard tool by Torrenc e & Webster! (Il999l) . We used the com- 
plex Morlet mother function, which is also suitable for estimat- 
ing the phase between two signals. In addition, we have also 
made use of the coherence spectrum. While the wavelet spec- 
trum gives information about the distribution of power in the 
frequency-time domain, the coherence of two time series, being 
dependent on the cross-wavelet transform, is useful for inves- 
tigating the interaction bet ween the two physical mechanisms 
JTorrence & Websterll 19991) : in our case between v/, and v.. 
The third and final step of our analysis represents an extension 
of the second step to the full population of magnetic features. In 
particular, we focused on the phase angle between v/, and v,. 
The total number of features tracked in the IMaX data and the 
simulation amounts to 2384 and 6388, respectively. The larger 
amount of features tracked in the simulation can be explained in 
terms of a combination of effects. Thus, the estimated average 
lifetime of the magnetic elements in the simulation is 2.5 - 3 
times smaller than the average lifetime of the magnetic features 
collected in the IMaX data. This fact, together with the larger 
duration of the simulation (58 min) with respect to the IMaX 
data (32 min), and the intrinsic higher spatial resolution of the 
simulation can significantly increase the number of features col- 
lected. 

Our goal is to obtain the distribution of the phase difference be- 
tween the transversal and longitudinal velocity inside the mag- 



netic elements. A clear phase relation would point to an inter- 
action between horizontal and vertical velocity perturbations. A 
number of constraints, described below, were applied on the fea- 
ture to ensure the reliability of the results. 
In the coherence diagram given by the wavelet analysis, we se- 
lected, for each magnetic feature, the location in the frequency- 
time domain at which the two time series under investigation 
show the highest coherence. The phase between the two signals 
is then estimated at this location. A phase estimate was consid- 
ered only if the coherence between the two signals under in- 
vestigation was above 0.8 outside the cone-of-influence of the 
wavelet phase diagram, i.e. in the region where the wavelet anal- 
ysis can be tiTisted. The temporal cadence of the simulations was 
reduced to match the IMaX cadence (33 s). 
Besides determining the phase shift between the longitudinal ve- 
locity Vj. and the transversal velocity v/,, we repeated the analysis 
between the surrounding oscillatory field v^nm and v/,. The latter 
is done to check whether any phase relation between the longi- 
tudinal and transversal velocity is strictly inherent to the wave 
signal observed inside the magnetic elements, or not. 
In addition, we also studied the distribution of periods at which 
V- and v/, are coherent. This is done to find the frequency band at 
which their interaction (if any) takes place. 
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Fig. 4. (a) Wavelet amplitude spectrum of V; from the simulation, (b) Wavelet amplitude spectrum of v;, from the simulation, (c) Wavelet phase 
coherence spectrum between v. and vi, from the simulation. The arrows represent the phase lag in polar coordinates (zero phase corresponds to 
arrows pointing to the right), (d) Wavelet amplitude spectrum of v^ from IMaX. (e) Wavelet amplitude spectrum of v;, from IMaX. (f) Wavelet 
phase coherence spectrum between v- and v;, from IMaX. The arrows represent the phase lag in polar coordinates (zero phase coiTesponds to 
arrows pointing to the right. The dashed line represents the cone-of-influence where results cannot be trusted (area between the dashed lines and 
the edges of the diagram). The white contours represent the 95% confidence level. 



4. Results 

4.1. FFT power spectra and amplitude of oscillations 

Through the features tracking we obtained the horizontal and 
vertical velocity associated to each magnetic element. 
As mentioned in Sect. 3.2, we at first estimated the average 
power spectra of the twenty magnetic features characterized with 
the longest lifetimes. 

In panels (a) and (b) of Fig. [3] we show the power spectra ob- 
tained from the IMaX data and the simulation. The power spec- 
trum for both horizontal velocity (panel a) and vertical veloc- 
ity (panel b) are shown. The vertical line indicates, the poorer 
frequency resolution among the time series chosen (i.e. the fre- 
quency resolution of the shortest time series v„„„ = 2/7'/f„^,/,). 
The region at low frequency below this line should therefore not 
be trusted. 

Both simulation and observations display a similar behaviour. 
The power spectrum of v~ shows that most of the power is con- 
centrated at low frequency, in the 2-7 mHz band. 
The horizontal velocity, in turn, is characterized by a broader 
power spectrum with high-frequency components (> 8-10 
mHz), although the highest power is still located in the low 
frequency band. The power of the horizontal velocity is larger 
than the power of v,, at all frequencies, particularly in the IMaX 
data. Consequently the amplitude of the v/, oscillations is larger. 
This can be seen in panels (c) and (d) of Fig. [3] where we plot 
the distributions of the standard deviation of v/, and v, obtained 
from the complete population of magnetic features. Both simu- 
lations and observations show an excellent agreement, with most 
of the magnetic features displaying horizontal oscillations hav- 



ing a mean amplitude around 1 - 1 .5 km/s, while for the longitu- 
dinal velocity, the corresponding values lie in the range 0.4 - 0.5 
km/s. 



4.2. Wavelet amplitude analysis and coherence: a case study 

Fig. |4] depicts the wavelet diagrams for the two longest time 
series found respectively in the IMaX data and in the simulation 
domain (i.e. the longest lived magnetic element in each). The 
time series associated with the magnetic element found in 
the IMaX data has been tiTincated to match the lifetime of 
the magnetic element from the simulation. This is to force 
the wavelet diagrams to be comparable in terms of frequency 
resolution. 

Panels (a) and (d) show the amplitude spectrum of Vj, for the 
simulation and the observations, respectively. Panels (b) and (e) 
show the amplitude spectrum of v/, and, panels (c) and (f), the 
coherence between v, and v/,. The arrows represent the phase 
angle between the two signals in polar coordinates, where zero 
is along the x - axis. In our sign convention, a negative phase 
angle (arrows pointing down) means that v/, leads Vj.. In the same 
plots, the white contour represents the 90% confidence level, 
while the dashed line represents the cone of influence. 
The vertical velocity, v-, is characterized in both the simulation 
and the IMaX data by oscillations with periods larger than 
130 - 150 s. It is worth noting that, although the fiducial 
line extends to longer periods in the three-minute band, only 
the region outside the cone-of-influence can be trusted and, 
therefore, only this region will be considered in our analysis. 
The horizontal velocity v/, of the magnetic elements show 
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Fig. 5. Power spectra of three different simulated random walks. The 
continuous hne represents the 95% confidence level, the dashed line 
represents the cone of influence. The Nyquist period is set by the sam- 
pling to 66 s. 



power in the band P = 100 - 150 s and at higher frequency. 
In particular in the simulation, where the Nyquist frequency is 
larger, one finds many peaks above the confidence level close to 
f = 50 s or even smaller periods. This result is compliant with 
what has found in the power spectral density (Sect. 4.1). 
In both the simulation and the IMaX data, one finds locations 
in the wavelet diagram where v, and v/, are highly coherent (see 
panels (c) and (f) of Fig. |4|i. As the coherence is a measure 
of the interaction between two signals, which depends on the 
cross-spectrum, one can find high coherence even outside the 
confidence levels of each time series. This means that even 
in the case of minimal power the processes associated to the 
vertical and horizontal velocity may still strongly interact. 
A high coherence is found to reside at periods of 100 s and 
shorter. The associated phase angle is always different from 
0°. In the case of the simulation the phase lag is in the range 
(-70°, -1 10°), while for IMaX data it is close to -95°. 



4.3. Power spectrum of random walks 

Small flux tubes in the photosphere are subject to the forcing ac- 
tion of the granular flows. Their displacements and oscillations 



represent therefore the response to this forcing. 
In this section we present the results of a simulation conducted 
to test whether the v/, power spectra are compatible with random 
walks. 

We simulated the paths of 500 2D random walks lasting for 2000 
s and with a temporal step of 1 s. The physical size of the spatial 
grid was set to 0.5 km. This is equivalent to set the amplitude 
of the fluctuation of the position in the time interval. This value 
was chosen a posteriori to ensure that the mean amplitude of 
V;, of the simulated random walks was ^1.3 km/s, thus similar 
to the real one at the same IMaX cadence (after resampling the 
simulated time series at 33 s). After setting up this value, we re- 
peated the simulation. For each random walk we estimated the 
power spectrum of v/,. In Fig. [3] (panel a) we plot (dot-dashed 
line) the average power spectrum obtained from averaging over 
all the 500 cases. The shape of the power spectra are similar, thus 
the observed oscillations of the flux tubes appear to be consistent 
with random walk displacements due to buffeting. For compari- 
son with the real case discussed in the previous section. In Fig. 
|5] we also show three examples of wavelet diagrams obtained 
from the simulated random walks. The power is distributed over 
a broad range of periods, from large periods at the frequency 
resolution imposed by the duration of the simulated time series, 
down to the Nyquist limit (66 s). 

4.4. Statistical analysis: histograms of phase 

To find a statistically reliable estimate of the phase relation 
^(Vj-v/,), we used the full sample of tracked magnetic elements. 
For this purpose, we used the wavelet analysis described in Sect. 
l4.2l to estimate the phase corresponding to the highest coherency 
in the wavelet diagram, for each tracked magnetic element. In 
Fig.|6^ we compare the results obtained from IMaX observations 
(continuous line) to those from the MHD simulation, resampled 
to match the IMaX cadence (dashed line). Both histograms dis- 
play clear maxima at -90° and H-90°. 

To test the consistency of this result, we studied the phase re- 
lation between the horizontal velocity and the longitudinal ve- 
locity surrounding the magnetic elements. The aim of this test 
is to check whether the ±90° phase angle is inherent to the ve- 
locity signal strictly inside the magnetic elements, or rather it is 
aff'ected by either contaminations due to the surrounding acous- 
tic field or some issue concerning the analysis. To this end we 
estimated the mean velocity signal in an annular region like that 
shown in Fig.[TJl, where we selected only the non-magnetic pix- 
els. The results of this check are plotted in Fig. |7] The phase lag 
is zero in this case. This result demonstrates that the +90° phase 
lag between the vertical and horizontal velocities is obtained 
only if the internal velocity field of the flux tube is considered. 
Therefore it is closely related to the nature of the MHD waves 
in the magnetic elements. As mentioned in the previous section, 
we tracked groups of pixels lying on the same 'hill' . This method 
tracks each magnetic element (peak in the Stokes V signal) sep- 
arately, even if it is part of a tight cluster (i.e. the magnetic signal 
stays high between individual elements. Clusters of magnetic el- 
ements may have a correlated horizontal velocity, hence their 
contribution to the histogram of phase may significantly affect 
its shape. To further test our results, we studied the histogram 
of phase obtained from the IMaX data but tracking contiguous 
pixels instead of local maxima. Thus, we considered all the con- 
tiguous pixels with Stokes V above a given threshold to belong 
to a single magnetic feature. In this case, regardless the presence 
of local maxima, the magnetic structure is tracked as a whole. 
The resulting histogram is also shown in Fig.|6^ (crosses). Even 
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in this case, the histogram of phase shows two maxima at -90° 
and +90°. 

We also note that a significant fraction of elements show a small 
phase shift (~ 23% in the IMaX case). This can be due to pro- 
jection effects. If a magnetic feature is inclined to the LOS, there 
will be a component of the transverse velocity, associated to 
the observed kink oscillations, along the LOS which will be: 
Vz.projected - Vhsin(y), whcrc y is the inclination angle. By the 
same token, the amplitude of the kink velocity will be decreased 
by cosy. 

If we consider the average observed horizontal velocity (see 
Fig. [3]l, after the correction for the projection effects, we get 
Vz,projected ~ 0.19 - 0.29 km/s. This is estimated using an in- 
clination angle between 10° and 15°, as obtained by Jafarzadeh 
et al. (in prep.) from IMaX observations. 

Form the histogram of the LOS velocity we can estimate the 
probability of finding a magnetic element with v/, smaller than 
the above limits. This corresponds to the case in which the lon- 
gitudinal velocity would be dominated by the projection of the 
horizontal velocity along the LOS, thus they would be in phase. 
We estimate this probability between 14% and 29%, which is 
comparable with the observed fraction of elements with a small 
phase shift. It has to be noted that our selection criteria based 
upon Stokes V and B^ in the IMaX data and the simulation, re- 
spectively, may result in a sample of magnetic features not re- 
stricted to vertical flux tubes, including a significant fraction of 
inclined features for which 0(v, - v/,) is close to zero because of 
the mechanism above. However, this may not be the only mech- 
anism responsible for the large spread of the peaks in the his- 
togram of phase. Another important cause may be the difference 
in the height of maximum response of Stokes V to vertical and 
horizontal velocities in the magnetic elements. This mechanism 
may contribute, along with the first one cited above, to increase 
the number of elements with small phase lag. 
Moreover, as anticipated in Sect. 3.1, we checked the results 
against possible effects due to the particular choice of the thresh- 
old used in the tracking code. In fig.|6]we also plot the results of 
this further analysis using a threshold of 2cr, 3cr and 4cr on the 
IMaX data. As expected the histograms are shifted to lower val- 
ues as the threshold increases, but their shape is not changed, 
showing again two maxima at -90° and -1-90°. 

We also studied the distribution of periods at which one 
finds the highest coherence between v, and v/,, that is periods 
at which the +90° phase lag is found. We recall here that we 
only considered those cases in which the coherence was at least 
0.8. This is a demanding constraint. This distribution is shown 
in Fig. |6}3. The continuous line represents the IMaX data, the 
dashed line the simulation. It is evident that, in most of the 
cases, the interaction between the vertical and the horizontal 
velocity takes place in the high-frequency band, at periods 
smaller than 200 s. 
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Fig. 6. (a) Histogram of phase ^(v- - v;,) obtained from IMaX with dif- 
ferent thresholds (2cr black continuous line, 3cr red dash-dot line and 
4(T dotted line), and the simulation (dashed line). Crosses represent the 
histogram of phase obtained from IMaX data tracking of contiguous 
pixels instead of local maxima, (b) Normalized histogram of the pe- 
riods at which the coherence, between V; and v/,, is largest for IMaX 
(continuous line) and the simulation (dashed line). 
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5. Discussion 



Fig. 7. Histogram of phase 0(vj - iy^m) obtained from IMaX data. 



The presented results concern the main properties of velocity 
perturbations in small magnetic elements in the solar photo- 
sphere. In particular we have studied transverse oscillations and 
longitudinal velocity oscillations of these magnetic elements. 
Starting from the longest time series associated with different 
magnetic features, we have retrieved the average power spectral 
density for the two kinds of oscillations. Since most of the 
magnetic features tracked, and their associated time series, 
have a short lifetime (shorter than 10 minutes), the wavelet 



analysis constitutes a more suitable tool for the analysis of the 
periodicities and, more specifically, for the estimate of the phase 
angle between v, and v/,. This is even more important when 
dealing with non-stationary signals like those expected for these 
kinds of velocity oscillations. 

Our results demonstrate that, while v, is characterized by 
comparatively low frequency oscillations (v < 7 - 8 mHz, i.e. 
periods larger than 120 - 140 s) and amplitudes of the order of 
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500 m/s, transverse displacements have a more extended power 
spectrum with very high frequency peaks (v ~ 8 - 10 mHz, i.e. 
periods of 120 - 100 s), and with amplitudes of the order of 1 
km/s. Both these modes are expected to be propagating at these 
frequencies with significant observed power. 
The wavelet analysis has also revealed that the horizontal 
velocity is characterized by high frequency wave packets, 
confirming the non-stationarity of the signals. Unfortunately, the 
short lengths of individual time series, due to the short lifetimes 
of the magnetic features, hampered the estimate of the mean 
duration of these packets, as well as the resolution of the low 
frequency part of the spectra. 

We have also demonstrated that at least the v/, power spectra 
associated to the longest lived magnetic elements are comp atible 
with those of random walks. iManso Sainz et al] (1201 Ih have 
shown that once small magnetic elements reach the intergranular 
lanes they are subject to the buffeting action at the hands of the 
surrounding granular and integranular turbulent flows and they 
follow random walks. However, a lot of photospheric magnetic 
elements may have ballistic trajectories before entering the 
intergranular lanes. Additionally, we used the coherence to 
investigate the interaction between the transverse velocity and 
the longitudinal perturbations. From the analysis we have found 
that in most of the cases a high coherence is found in the high 
frequency band (P < 200 s). This demonstrates that, at high 
frequency, the vertical and the horizontal velocity are strongly 
coupled and interacting. This interaction is also characterized by 
a phase lag. The statistical analysis has shown, in fact, that the 
phase difference 0(v, - v/,) is found to be ±90°. It is not possible 
to solve the ambiguity which is inherent to our methodology. 
Note that only those cases in which the coherence was above a 
very high threshold (0.8) were considered, thus strengthening 
the reliability of the results. 

This result demonstrates that buffeting-induced displacements 
of the flux tubes are accompanied by, and possibly excite, 
longitudinal and compressive MHD waves. Jafarzadeh et al. 
(in prep.) have found, using multi-wavelength observations, 
obtained from Sunrise/SuFI, that small magnetic elements are 
mainly dominated by fast MHD waves with velocity of the order 
of 40 km/s, as well as kink waves with velocities in the same 
range. Both kinds of waves were upward propa gating. This 
scena rio is also supported by numeri c al simu lations (iNutto et alJ 
l2012h . On the other hand. lJess et al.l (1201 2bl) have also detected 
slow upward propagating longitudinal waves in small magnetic 
elements, using high cadence broad band 2D data, visible as 
periodic intensity fluctuat i ons in the range 1 10 - 600 s. 
iNakariakov & Verwichtel (l2005h have reviewed in detail the 
excitation of waves in a magnetic cylinder. They have shown, 
by means of numerical simulations, that kink waves and 
longitudinal compressive waves may coexist in the same flux 
tube, the latter being excited by the horizontal motion of the 
flux tube itself (see for example their Fig. 4 and the associated 
movio- In particular, the kink mode excites a longitudinal 
density perturbation, whose maximum is reached at the inver- 
sion points of the flux tube's horizontal motion (i.e. when the 
horizontal velocity is zero). The longitudinal perturbations and 
the horizontal velocity are therefore 90° out of phase. Although 
our results are consistent with this scenario, we believe that the 
phase shift alone is not enough to uniquely identify the mode 
of oscillation. This is because, the vertical velocity in both 
the observations and the realistic MHD simulations used, is 



probably the result of a mixture of modes, whose superposition 
hinders their identification. 

The interaction and, therefore, the energy exchange between the 
kink mode and the longitudinal mode takes place over a wide 
range of frequencies (see Fig.|6]panel b), being more probable at 
high frequencies {P < 100 s). At the frequencies of interaction, 
the power spectra of both kink oscillations and longitudinal 
oscillations show less power than the low-frequency band. 
It is worth stressing that, transverse perturbations and the 
longitudinal mode can interact even where the power is below 
the maximum power achieved in the power spectrum associated 
with each time series. In other words, in order for the two modes 
to interact, it is not necessary to have a large amount of power 
at the same time and at the same frequency in both processes. 
Although the amount of energy per unit frequency is small in 
this case, Fig.|6]demonstrates that this interaction can take place 
over a wide range of frequencies (possibly even at frequencies 
above the Nyquist frequency imposed by the cadence of our 
observations), thus the total amount of energy involved could be 
significantly larger. 

These results strongly rely on the coherence and the phase 
estimates. Unfortunately, while for the power it is straightfor- 
ward to give an estimate of the confidence level of its peaks in 
the wavelet diagram, this is not true for the phase. The major 
problem comes from the fact that, in the case of pure noise, there 
is no reference value for the angle, the latter being a variable 
uniformly distributed from -n to +7t. For this reason, we have 
provided a statistical check of our phase estimate, by exploiting 
the large number of magnetic elements at our disposal. This not 
only proves the reliability of the phase relation between V; and 
v/,, but also underlines the general character of this result, which 
holds for a significant fraction of magnetic elements. 
Moreover, we have also checked whether this phase relation 
is found only when considering the velocity field inside the 
magnetic features or, rather, it is obtained also when the vertical 
velocity outside the magnetic elements (vlnm) is considered. 
This check has demonstrated that the ±90° phase can be only 
obtained in the first case. This basically means two important 
things. Firstly, our phase estimate is not affected by contam- 
ination through the acoustic field surrounding the magnetic 
elements. Secondly, the results seen so far are closely related 
to the MHD waves present in small magnetic features and are 
not significantly affected by spurious effects coming from the 
external velocity field. This means that we are really taking 
advantage of the high spatial resolution provided by SUNRISE, 
which allows us to isolate the oscillation field inside small-scale 
magnetic elements. 

We compared the observational results with radiative-MHD 
simulations, and we found them in excellent agreement. This 
demonstrates the reliability of the IMaX results and puts our 
findings on a firmer ground. 

The phase histograms could suffer from the presence of a bias 
due to the inability to track magnetic features for a long time. 
This is most likely due to the short lifetimes of the magnetic 
features, although we cannot rule out detection problems in 
some cases (e.g. when a magnetic feature is lost or confused 
due to interactions with other features). Therefore, longer lived 
magnetic features could host a lot more power in longer-period 
oscillations. From the theoretical point of view, the presence 
of detectable longitudinal velocity oscillations inside the flux 
tubes, associated with buffeting-induced transver se perturba- 
tions, is in agreeme nt with the theoreti c al work by iHasan et al] 
(2003) and Musielak & Ulmschneider (2003b). These authors 



' The movie of the simulation is available online at the following link: 
[http : //solarphysics . livingreviews . org/Articles/lrsp- 2(S)(S)5-ty/q^i««&gc^cyMriH^a8B^/S^^$a^%Mng}^/;W^atJj<^ 
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due to the external forcing by granular buffeting, not only 
generates transverse oscillations, but also longitudinal oscil- 
lations, which should then be observable inside the magnetic 
structures as Doppler velocity oscillations. This scenario is 
in good agreement with our findings. Moreover, our results 
give new insight into the interaction between longitudinal and 
transverse perturbations and, as far as we know, constitutes the 
first observational evidence of the interaction between these two 
in solar magnetic features. To this regard, one way to distinguish 
random horizontal motions of flux tubes from (propagating) 
kink waves would be to look at different heights simultaneously. 
In this case, for not vanishing phase lags this would indicate the 
presence of propagating waves. One may think that if a tube is 
displaced by convection in the deeper layers in one horizontal 
direction, it will move first in the lower layers then in the upper 
layers, mimicking a propagating wave. However, such a motion 
may set up a propagating kink wave. 



6. Conclusions 

In this paper we have reported the main properties of MHD 
waves and perturbations in small-scale magnetic elements by 
comparing high-resolution SUNRISE/IMaX observations to 
MHD simulations. This analysis reveals the interaction between 
transverse perturbations and longitudinal velocity oscillations in 
small-scale magnetic elements. 

We find small-scale magnetic features in the solar photosphere 
to host a rich variety of MHD perturbations which we studied by 
means of a statistical approach. Our analysis has demonstrated 
that most of the magnetic elements observed display transversal 
oscillations with amplitudes of the order of 1 km/s, and with a 
power spectrum which extends to very high frequencies (v > 10 
mHz). The longitudinal waves found in the data are character- 
ized by a lower frequency (v < 7-8 mHz). Although their 
power spectra present peaks at different positions, the interac- 
tion between kink displacements and longitudinal waves is, with 
high confidence level, found to take place in the range P < 200 
s. In particular we found that the Doppler velocity signal within 
the magnetic features shows a preference for +90° phase shift 
with respect to the horizontal velocity of the magnetic feature 
itself. 

This result is only obtained if the velocity signal inside the mag- 
netic elements is considered, thus ruling out the contamination 
from the external non-magnetic oscillatory field. 
These results are compliant with the theoretical framework in 
which kink waves, generated by the granular buffeting, are ac- 
companied by longitudinal waves which are generated through 
non-linear interactions. 

Small scale magnetic fields cover a significant fraction of the so- 
lar surface, therefore the presence of this kind of MHD waves 
and their energy exchange, may assume a remarkable impor- 
tance in the energy transportation toward the upper layers of 
the Sun's atmosphere. This work has profited greatly from the 
high resolution and the absence of seeing characterizing the 
SUNRISE data, which make them ideal for studying time- 
dependent phenomena. New high spatial resolution and simul- 
taneous multi-line spectropolarimetric observations of the solar 
quiet photosphere and chromosphere are needed to reveal clues 
about the propagation of these MHD waves in small scale flux 
tubes. 

Furthermore, our results emphasize the need for very high tem- 
poral cadence (< 30 s) data in order to resolve the very high 
frequency part of the power spectrum. 
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